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ABSTRACT 

We suggest a solution to an important problem of observational helioseismology of 
the separation of lines of solar acoustic (p) modes of low angular degree in oscillation 
power spectra by constructing optimal masks for Doppler images of the Sun. Accurate 
measurements of oscillation frequencies of low-degree modes are essential for the deter- 
mination of the structure and rotation of the solar core. However, these measurements 
for a particular mode are often affected by leakage of other p modes arising when the 
Doppler images are projected on to spherical-harmonics masks. The leakage results in 
overlaping peaks corresponding to different oscillation modes in the power spectra. In 
this paper we present a method for calculating optimal masks for a given (target) mode 
by minimizing the signals of other modes appearing in its vicinity. We apply this method 
to time series of 2 years obtained from Michelson Doppler Imager (MDI) instrument on 
board SOHO space mission and demonstrate its ability to reduce efficiently the mode 
leakage. 

Subject headings: Sumoscillations 
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1. Introduction 

Accurate measurements of frequencies and other 
physical properties of solar oscillations are crucial for 
understanding the internal structure and dynamics of 
the Sun. Up to now the determination of low-degree 
p-mode parameters like cyclic frequency v, linewidth 
or frequency splitting is based on fitting the so-called 
m-v diagrams obtained by projecting time series of 
solar images (usually Dopplcrgrams) on to a basis of 
spherical harmonics of angular degree / and azimuthal 
order m. Because only one hemisphere of the Sun 
is observed the projection on to spherical harmonics 
is not orthogonal. This results in leakage of modes 
which makes the line fitting quite complicated espe- 
cially when modes overlap with the target mode. We 
have therefore decided to apply a method originally 
developed by Kosovichev (1986) based on a singu- 
lar value decomposition (SVD) technique in order to 
minimize the leakage. 

We discuss in Sec. 2 the principle of the method, 
and introduce the concepts of 'global' and 'local' op- 
timal masks. The global optimal masks are designed 
to minimize the contribution of all modes except the 
target one. The local masks minimize only modes 
within a narrow frequency interval centered at the tar- 
get mode. In Sec. 3 we show the theoretical efficiency 
of the method compared to the usual spherical har- 
monics projection. In Sec. 4 we illustrate this method 
with a 2-year time series of MDI LOI-proxy data, ob- 
tained by rebinning of the original 1024 by 1024 CCD 
pixels into 180 bins (Hoeksema et al., 1998; Scherrer 
et al., 1995) according to the shape of the LOI instru- 
ment (Appourchaux et al. 1997). 

2. The optimal-mask method 

The aim of the optimal masks (Gough, 1984 and 
Christensen-Dalsgaard, 1984) is to find an optimal 
linear combination of N signals arising from N bins of 
the solar image in order to produce a time series con- 
taining mainly the contribution of a mode with given I 
and m, the so-called target mode, and minimizing the 
contribution of others modes. It is particularly im- 
portant to minimize the contribution of modes which 
are close to the target mode in the power spectrum. 
These are mainly modes of the same rotationally split 
multiplct as the target mode. At frequencies about 
2.5 mHz and higher the linewidth of p modes becomes 
larger than the separation between the peaks. There- 
fore, it is important to separate the overlapping modal 




Fig. 1. — An example of MDI LOI-proxy Doppler- 
grams. The full-disk 1024xl024-pixel MDI Dopplcr- 
grams are rebinned on board SOHO into 180 bins 
with the boundaries along heliosphcric meridians and 
parallels. The variation from dark to light colors is 
due to solar rotation. 

lines. We describe this method in the application to 
the MDI LOI-proxy data. However, this method is 
quite general and can be applied to other spatially 
resolved helioseismic data. 

Because we observe only one hemisphere of the 
Sun, the commonly used spherical-harmonic masks 
are not able to remove all the components of the tar- 
get mode multiplet. Instead of using the spherical 
harmonics masks we apply a SVD method to derive 
the optimal masks. To apply this method we first 
need to model the signal produced in each image bin 
by a mode of given angular degree I, azimuthal order 
m and radial order n. For simplicity, in this paper, 
the horizontal component of velocity is neglected be- 
cause we are only interested in low-degree p modes. 
The inclusion of the horizontal component is straight- 
forward, and generally improves the quality of the 
optimal masks in the low-frequency part of the solar 
oscillation spectrum where this component becomes 
significant. Hence, to the first order, the signal of a 
normal mode in each pixel j of the CCD detector is: 

v { JL(t)=^{v nlm e^s<lil(t)}, (1) 
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where Vnim is the mode amplitude, u> n im is the mode 
frequency, and s^t) = Y im [0* (t) , <j)* (t)]nj is the sen- 
sitivity function which depends on the spatial eigen- 
function, spherical harmonic Yi m , and the factor jij 
arising from the projection of the radial velocity onto 
the line of sight. Both Yj m and fij are expressed in 
terms of the heliographic coordinates, 0*(t) and 4>*j(t), 
of the pixels, which are generally functions of time t. 
We have to take into account the fact that the orien- 
tation of the solar image, which is usually described 
by two angles of the rotation axes with respect to 
the ecliptics, Po and Bq, and the image size on the 
detector are possibly changing with time. The time 
dependences of the P and P angles should there- 
fore be taken into account by rotating the spherical 
coordinates accordingly. In our case the variations of 
the image size were not taken into account because 
the data were combined into fixed 180 bins on board 
spacecraft prior to any mask application. An example 
of the rebinned Dopplergrams ('MDI LOI-proxy') is 
shown in Fig. 1. 

The k th signal of each of the N bins is obtained 
by averaging the signals coming from Nk CCD pixels 
withing this bin as defined by the mapping between 
the MDI pixels and the LOI-proxy bins: 
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Ideally, if the sum in the square brackets is equal 
to the product of two Kroneker symbols, Su S mmo , 
then v no i omo (t) consists only of the signal of the tar- 
get mode. However, in practice, this sum does not 
represent 5u a 5 mma and contains not only the target 
mode but also contributions from other modes. 
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optimal mask if this vector is chosen to maximize the 
signal for the target mode and minimizes the signals 
for other, say M, modes located close to the target 



mode in the Fourier domain. These conditions can be 
written as a minimization of the quadratic form: 

r = ||S.w-I|| 2 , (4) 

where S is a (M + 1) x N matrix with elements 
S^i m {t) and I is a (M + l)-element vector with el- 
ements Su 5 mmo , l a and m being the angular degree 
and order of the target mode. 

Finding w is equivalent to finding a least-squares 
solution to the system of linear equations, S • w = I. 
Depending on M, the number of modes we want to 
minimize, we have a system which is over-determined 
(N > M+l) or under-determined (N < M+l). For 
both cases the problem is efficiently solved using the 
singular value decomposition (SVD) of the matrix S 
(e.g. Press et al. 1992). This method optimally solves 
the linear system and its solution w gives the optimal 
mask. This mask provides an efficient filtering of the 
selected M modes. Of course, the minimization is less 
efficient when M is large. There is therefore a trade- 
off between the number of modes to minimize and the 
efficiency of this minimization. 

Matrix S describes the sensitivity of each detec- 
tor bin to the mode of oscillation. For the special 
case of an under-determined system (N < M+l) this 
technique is equivalent to the technique of diagonal- 
ization of the leakage matrix (Appourchaux et al., 
1998). The SVD method allows us to consider gen- 
eral cases. With this technique we can specifically 
decide which modes must be filtered and to what ex- 
tend their contribution to the power spectra should 
be reduced, making the technique extremely flexible. 
In practice we filter as much as possible the modes 
having frequencies falling in a given frequency win- 
dow around the target making a kind of local cleaning 
of the spectra. We therefore call the corresponding 
masks 'local optimal masks' to distinguish from the 
'global optimal masks' obtained when minimizing the 
contribution of all the modes not only those around 
the target mode. The larger the window around the 
target the more modes have to be minimized, and thus 
the less efficient is the filtering. In principle we are 
able to remove efficiently as many modes as the num- 
ber of the binned signals, 180 in our case. Above that 
number the SVD acts as a least square minimization 
and the modes cannot be fully filtered. In practice, 
the efficiency of the filtering depends on the selected 
mode set and the noise level. 

The solar noise obviously will be less filtered than 
modes by the masks because its spatial structure is 
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incoherent. Thus the signal-to- noise ratio of the tar- 
get mode depends on the mask. More efficient masks 
typically require higher weights w^ m , and, therefore 
amplify the noise. This leads to the trade-off between 
the mask efficiency and the noise level. To adjust this 
trade-off we slightly modify Eq. (4) and add the noise 
contribution through a regularization term. Then the 
quadratic form to minimize is 

r = |S-w-I|| 2 + a||n-w|| 2 , (5) 

where n = {n^\ nW} is the vector of the stan- 
dard deviation of the noise in each of the N bins, and 
a is the regularization parameter. The SVD tech- 
nique is still applicable here. Eq. (4) corresponds to 
the case of a=0 which means no regularization. 

The noise level is estimated by taking the standard 
deviation of the signal in each bin, assuming the main 
contributor is the solar noise. Generally, the noise 
level is the lowest at the center of the solar disk and 
increases toward the limb. Table f gives a sample of 
the relative noise level of the MDI measurements as a 
function of the angular distance from the disk center. 
Obviously, the absolute noise level is not required for 
Eq.(5). 

The regularization parameter a is chosen such as 
the amplitudes of the filtered modes is reduced to 
the level of noise, ft turns out that the signal-to- 
noise ratio obtained with the local optimal masks is 
similar to what is obtained with spherical-harmonic 
masks, though it may be somewhat smaller (up to 
50%). Only for /=0 the difference in the signal-to- 
noise ratio between the two masks can reach a factor 
of two or more. Because of solar noise we also avoid 
using the limb bins. These bins improve very little the 
efficiency of the masks but contribute significantly to 
the noise. Removing not only limb bins but also the 
next round of bins helps to increase the signal-to- noise 
ratio at low frequencies for I = modes, but this is 
not effective for modes of higher degrees. 

The power of this technique is also that the optimal 
masks for any P and B angles are easily obtained 

Table 1 : Relative noise level, n, in the MDI data as 
a function of the angular distance, d, from the disk 
center. 



d, deg. 1.72 11.0 20.5 30.7 42.1 56.1 
n 1.00 1.15 1.40 1.73 1.99 2.27 



applying the spherical harmonic rotation matrix to 
the optimal masks with Po~0 and B o =0- Therefore it 
is possible as for MDI to apply the optimal masks for 
each day according to the value of Bq which is chang- 
ing with time without recomputing the SVD problem 
which is a heavy computational task . 

3. Comparison of spherical-harmonic masks 
and optimal masks 

In order to test the efficiency of the optimal masks 
we compute how these masks filter the modes ver- 
sus the target mode and compare to the spherical- 
harmonic masks. We assume that each mode leads to 
a velocity signal which is a perfect spherical harmonic 
projected onto the line of sight. The degree range of 
modes is given by the properties of the detector. The 
180-bin detector (the LOI-proxy) shown in Fig. 1 is 
sensitive up to modes 1=15. The higher-degree modes 
do not contribute significantly to the observed signal. 
Therefore, in the global masks we restrict the mode 
range from I = to 15. 

The results of the global and local masks are rep- 
resented in the form of a map (called hereafter am- 
plitude map) showing the relative amplitude of the 
modes identified with the doublet (I, m) for I between 
and 6 and m between —I and +1. 

Figure 2 shows, as an example for a target 1=1 
and m=l, the way various masks filter the unwanted 
modes. On top left we used the disk-integrated signal 
which leads to the usual problem of mode blending 
between m=± 1. It is well-known (Chang, 1997) that 
this blending especially bias the determination of fre- 
quency splitting. For the spherical-harmonic masks 
(top right), the blending is smaller but still there is 
the same problem. Using the local optimal masks for 
a frequency window of ±15 /iHz we are able to reduce 
below the noise level the m=-l as well as the /=3,6 
and 9 modes which were also interfering with the tar- 
get mode. Note that global masks (bottom, right) 
as expected do not remove perfectly these modes be- 
cause the number of modes to filter (around 3000) 
exceeds the number of bins by a large factor. We 
give another example for a l=A, m=3 mode in Fig. 
2. Applying the local optimal masks should therefore 
lead to clean power spectra around the target mode 
as shown in next section using the MDI data. 
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Fig. 2. — Amplitude maps for a disk-integrated 
mask (top, left), spherical- harmonics masks(top, 
right), local optimal masks (bottom, left) and global 
optimal masks (bottom right) with target 1=1, m=l. 
The size of the symbols on this map is proportional 
to mode amplitudes. 
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Fig. 4. — Power spectrum for 1=1 m=-l mode using 
the spherical-harmonic mask (grey curve) and local 
optimal masks (black curve). 
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Fig. 3. — The same as in Fig. 2 but for 1=4, m=2. 
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Fig. 5. — Power spectrum around target mode 1=0 
at 3034 /xHz for local optimal masks with a different 
regularization parameter a: zero (black solid curve), 
small (dotted curve), large (dashed curve) and for in- 
tegrated velocity (grey curve). 
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4. Optimal masks applied to MDI data 

In order to apply the optimal masks to MDI data 
we need in principle to take account of the time- 
dependence of the signal due to the displacement of 
the image on the detector and change of its size as 
the SOHO-Sun distance is changing. However, for the 
low-degree modes this effect is expected to be small, 
and was neglected in these studies. 

The larger effect would result from the variations of 
the images orientation angles Po and Bq. However, in 
the MDI observations the Po is quite stable and much 
less than 1 degree. The Bo varies between -7 and + 7 
degrees during a year. However, we found that these 
variations lead only to non-significant variations of 
our optimal masks. Therefore, we adopt the mean 
value, B o =0. 

As an example, Figure 4 shows overlapped two 
pieces of power spectra of a 1=1, m = —1 mode, 
one obtained with the standard spherical harmonic 
masks (grey curve) and the other one with the lo- 
cal optimal masks (black curve). In the first case, 
the power spectrum shows a significant leakage of the 
1=1 m=l mode which belongs to the same rotation- 
ally split multiplct. This mode partially overlaps with 
the target mode. In the optimal mask spectrum, the 
contribution of this mode is reduced to the level of 
noise. As previously noticed the signal-to-noise ra- 
tio is better with the spherical-harmonic masks but 
the leakage mode is strongly blended with the target 
one leading to bias in the estimation of mode param- 
eters. Modes of Z=6 and 9 were also filtered for this 
case. The singular value cut-off was chosen to be of 
the order of the numerical accuracy and no regular- 
ization was applied which means optial filtering with 
a slightly larger noise level. 

Figure 5 shows the effect of regularization for a 
target mode of 1=0, where the modes of 1=2, 5, 8, 
11 and 15 are filtered out in a window of ±15yuHz 
around the target mode. For a=0 the contiguous 1=2 
mode multiplet at 3024 /xHz is reduced to noise level 
whereas as a increases this mode becomes larger and 
the noise level decreases. We found that in most cases 
the optimal masks with a = work sufficiently well. 

5. Conclusions 

Mode leakage which results in overlapping mode 
peaks in solar oscillation power spectra is one of the 
most significant problems in observational hclioseis- 



mology. It may lead to significant errors in the deter- 
mination of mode frequencies and other properties. 
We have shown that this problem can be efficiently 
solved by using the optimal masks which are applied 
to series of solar images and reduce the signals of so- 
lar modes in a narrow frequency interval around a 
target mode to the level of noise. The optimal mask, 
however, may increase the noise level. Therefore, the 
filtering efficiency of the masks has to be balanced 
with the noise level. This is achieved by a regular- 
ization technique. We have demonstrated in the case 
of low-degree modes that the optimal mask technique 
allows us efficiently isolate individual modes in rota- 
tionally split multiplets, which are unresolved in the 
commonly used disk-integrated data, and also isolate 
modes in the case of overlapped multiplets of different 
angular degree. The optimal masks are quite efficient 
for filtering out unwanted modes as long as their num- 
ber is less or of the order of the number of solar image 
bins. We used the MDI LOI-proxy data to calculate 
the optimal masks. However, this method is quite 
general and can be used for any spatially resolved he- 
lioseismic data. 

The method has been applied for the determina- 
tion of the central low-degree frequencies by Toutain 
et al (1999). In a future paper, we will present the 
results for rotational splitting of low-degree modes. 
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